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pH ■ Abstract 

We investigate the inclusive top quark pair production in proton-proton and proton- 
*> I antiproton collisions at high energies in the framework of fcr-factorization QCD approach. 

jj-s ' Our study is based on the off-shell partonic subprocesses g*g* — )■ QQ and q*q* — )■ QQ, where 

lO ■ the transverse momentum of both incoming quarks and gluons are taken into account. The 

unintegrated parton densities in a proton are determined by the Kimber-Martin-Ryskin 
prescription as well as the CCFM evolution equation. The conservative error analisys is 
performed and comparison with the results of traditional NLO pQCD calculations is done. 
Our predictions agree well with the recent experimental data taken by the DO, CDF, CMS 
and ATLAS collaborations at the Tevatron and LHC energies. 
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With the startup of the LHC, high energy particle physics entered a new era. One of 
the important part of the LHC program are the investigations of top-antitop pair produc- 
tion [HIS]. These studies are interesting for several reasons. At higher (compared to the 
Tevatron) energies, measurements with higher precision become available, which provide an 
opportunity to test the predictions of various theoretical models. The top quark mass rrit 
is an important parameter of the Standard Model (SM) and it affects predictions of other 
SM observables via radiative corrections. A precise measurement of rrit is crucial since it is 
one of the most important inputs to the global electroweak fits which provide constraints on 
the SM itself, including indirect limits on the mass of the undiscovered Higgs boson. More- 
over, top pair production is an important background in various searches for physics beyond 
the SM. Due to the large mass of the top quark, many models of physics beyond the SM 
predict observable effects in the top quark sector which can affect the top quark production 
rate. Finally, the abundant ti sample can be used for improving many aspects of detector 
performance. 

From the theoretical point of view, ti quark pairs in pp and pp collisions are produced 
via standard QCD gluon-gluon fusion and quark-antiquark annihilation subprocesses. At 
the LHC, the production mechanism is dominated by the gluon-gluon fusion, whereas at 
the Tevatron, top quark pairs are predominantly produced through quark-antiquark anni- 
hilation. In the framework of standard QCD, ti pair production at modern energies has 
been studied in many papers (see, for example, [3H9] and references therein). In particular, 
theoretical calculations of ti cross section have been carried out at next-to-leading order 
(NLO) [TIIH] and approximate next-to-next-to-leading order (NNLO) [9] of QCD. The re- 
sults of these calculations agree with the Tevatron and LHC data within the theoretical 
and experimental uncertainties. However, an alternative description can be provided by the 
/cT-factorization approach of QCD [lOJ. A detailed description and discussion of the kr- 
factorization formalism can be found, for example, in reviews [11]. Here we only mention 
that this approach has been successfully applied recently to describe the heavy flavour [T^1IT3] 
and prompt photon |14j production at the Tevatron and LHC energies. 

The present note is motivated by very recent measurements [Tl[2] performed by the CMS 
and ATLAS collaborations where the ti pair cross section in pp collisions at ^/s = 7 TeV 
has been measured. Our main goal is to give a systematic analysis of first LHC data [Tl[2] 
as well as recent Tevatron data [T514T7] in the framework of fc-r-factorization approach. We 
consider both gluon-gluon fusion and quark-antiquark annihilation subprocesses and take 
into account the non-zero transverse momentum of incoming gluons and quarks in a proper 
way. Specially we study different sources of theoretical uncertainties. Such calculations are 
performed for the first time. 

The cross section of ti pair production in the /cT-factorization approach is calculated as 
a convolution of off-shell (/cT-dependent) partonic cross sections and corresponding uninte- 
grated parton distributions in a proton. The contributions from the gluon-gluon fusion and 
quark-antiquark annihilation can be presented in the following simple forms: 
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where /g(x, k|.,/i^) and /^(x, k^,//^) are the unintegrated quark and gluon distributions in 
a proton, \M{g*g* -)• tt)p and |7W(g*g* -> tf)p are the off-shell (depending on the k^^^ 
and k22^) matrix elements squared, and s is the total center-of-mass energy. The produced 
top quark and anti-quark have the transverse momenta piT and p2T and the center-of- 
mass rapidities yi and 1/2. The initial off-shell partons have fractions xi and X2 of parent 
protons longitudinal momenta, non-zero transverse momenta ki^ and k2T and azimuthal 
angles 0i and 02- The evaluation of \M.{g*g* — )■ QQ)P, where Q is any heavy quark, is 
described in detail in [18]. The evaluation of \M.{q*q* — )■ QQ)\^ is straightforward. Here 
we only mention two technical points. First, according to the fc^-- factorization prescription 
|10j . the summation over the incoming off-shell gluon polarizations in (1) is carried with 
Y^ e'^e" = kj.k^/k|., where ky is the gluon transverse momentum. Second, when we calculate 
the squared matrix element in (2), the spin density matrix for off-shell spinors is taken in the 
form u{p)u{p) = xpp [19j, where x is the fraction of initial proton longitudinal momentum 
Pp. In all other respects our calculations follow the standard Feynman rules. Note that since 
the expression for the off-shell quark spin density matrix has been derived in the massless 
approximation, numerically we neglect the masses of any incoming quarks. 

To determine unintegrated quark and gluon densities in a proton we apply the Kimber- 
Martin-Ryskin (KMR) approximation [20] . This approach is the prescription to construct the 
unintegrated parton distributions from the known conventional ones. In this approximation, 
the unintegrated quark and gluon distributions are given by 
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where Pab{z) are the usual unregulated LO DGLAP splitting functions. The theta func- 
tions which appear in (3) and (4) imply the angular-ordering constraint A = yu/(/i + Ik^l) 
specifically to the last evolution step to regulate the soft gluon singularities. The Sudakov 
form factors Tg(k|,,/i^) and Tg(k|,,/i^) enable us to include logarithmic loop corrections to 
the calculated cross sections. Numerically, for the input we have used leading-order parton 
densities g(x,/i^) and g{x,fi'^) from recent MSTW set pT] . 

Another solution for the unintegrated gluon distribution have been obtained [22] from 
the CCFM evolution equation where all input parameters have been fitted to describe the 
proton structure function F2{x,Q'^). Below we will use proposed gluon density (namely, set 
AO) to study the uncertainties of our calculations connected with the non-coUinear evolution. 



Since we take into account the quark-antiquark annihilation subprocess, we should apply 
relevant unintegrated quark distributions to accomplish the CCFM gluon ones. At present, 
however, unintegrated quark densities are available in the framework of KMR approach onljo, 
therefore we will use following simple approximation. We will divide the quark densities into 
several parts which correspond to the interactions of valence quarks, sea quarks appearing 
at the last splitting and sea quarks coming from the earlier (second-to-last, third-to-last and 
other) steps of gluon evolution. First of them, unintegrated valence quark distributions have 
been obtained [23] from the numerical solution of the CCFM-like equatioiicl. To calculate the 
contribution of the sea quarks appearing at the last step of the gluon evolution we convolute 
the CCFM-evolved unintegrated gluon distribution /^(x, k|,,yU^) with the standard leading- 
order DGLAP splitting function Pqg{z): 



/f '^)(x,4,/i^) = ^^//,(x/^,k^,/i2)P,,(z)rfz. (5) 

X 

Note that in the region of small kf. < Qq the scale in the strong coupling constant a^ is kept 
to be fixed at go = 1 GeV. To estimate the contribution of the sea quarks coming from the 
earlier evolution steps we apply the procedure based on the specific properties of the KMR 
scheme. Modifying (3) in such a way |2l] that only the first term is kept and the second 
term is omitted and keeping only the sea quark in first term, we remove the valence and 
j{sea)i^^^ ky, jj?) quarks from the evolution ladder. In this way contribution of the sea quarks 
coming from the earlier gluon splittings to the /g(x, k^,/i^) is taken into account only. 

Other essential parameters were taken as follows: top quark mass mt = 170 GeV, renor- 
malization and factorization scales /z = ^rrit (where we vary the parameter ^ between 1/2 
and 2 about the default value ^ = 1 in order to estimate the scale uncertainties of our cal- 
culations), LO formula for the strong coupling constant ^^(/i^) with Uf = 4 massless quark 
flavours and Aqcd = 200 MeV, such that as{M'^) = 0.1232. The multidimensional integra- 
tion in (1) and (2) has been performed by means of the Monte Carlo technique, using the 
routine VEGAS ^25j. The full C-|--|- code is available from the authors on requescl. 

We now are in a position to present our numerical results. In Fig. 1 we confront the cal- 
culated top quark pair total and differential (as a function of t-quark transverse momentum) 
cross sections in pp and pp collisions with the recent data p!|[2l [T5] - [T7t l26 t [27] taken by the DO, 
CDF, CMS and ATLAS collaboration at the Tevatron and LHC. For comparison we plot 
also the NLO pQCD predictions ^ listed in [15]. One can see that the experimental data 
are reasonable well described by the /cy-factorization approach. Concerning the transverse 
momentum distribution, our predictions tend to overshoot the DO data (as well as the NLO 
pQCD predictions) at high pt values, but still agree with data within the scale uncertain- 
ties. The latter, since we are working in the LO approximation, are larger compared to the 
uncertainties of NLO pQCD calculations, which are of about 8 or 10% [ll[2]. The evaluation 
of higher-order corrections in the /cT-factorization approach is the special task and is out 
of our consideration. We would like to note, however, that the fcT-factorization approach 



^See, for example, reviews [TT] for more details. 
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distributions. 
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at LO level automatically incorporates the main part of the standard (collinear) high-order 
corrections yjj. The dependence of our predictions on the unintegrated parton densities 
is rather weak: the results coming from the CCFM and KMR parton densities practically 
coincide at the Tevatron and are very similar to each other at the LHC conditions. 

The relative contributions of gluon-gluon fusion and quark-antiquark annihilation sub- 
processes are shown in Fig. 2, where we used the KMR parton densities for illustration. 
Specially we estimate the contribution of the sea quarks coming from the not last gluon 
splittings (dotted curves) and find it to be of about 10%. Therefore, the uncertainties con- 
nected with the unintegrated parton densities are smaller than the scale uncertainties of our 
calculations. To avoid the possible double counting the contribution of the sea quarks coming 
from the earlier gluon evolution steps are not included in the CCFM predictions shown in 
Fig. 1 since part of them can be already included into the CCFM results (via initial parton 
distributions which enter to the CCFM equation). Note that these contributions are already 
included in the quark-antiquark annihilation predictions in Fig. 2. 

To conclude, in the present note we apply the CCFM and KMR unintegrated parton 
densities to the analysis of the first experimental data on the top quark pair production taken 
by the CMS and ATLAS collaborations at the LHC and the recent data taken by the DO 
and CDF collaborations at the Tevatron. Using the off-shell matrix elements of gluon-gluon 
fusion and quark-antiquark annihilation subprocesses (where the transverse momentum of 
both incoming quarks and gluons are properly taken into account) we have obtained a 
reasonably good agreement between our predictions and the data. It is important for further 
studies of small-x physics at hadron colliders, and, in particular, for searches of effects of 
new physics beyond the SM at the LHC. 
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Figure 1: The total and differential cross sections of the top quark pair production in pp and 
pp collisions at the Tevatron and LHC The solid and dash-dotted histograms correspond to 
the results obtained with the KMR and CCFM parton densities, respectively. The upper 
and lower dashed histograms correspond to scale variations in KMR predictions, as it is 
described in the text. The dashed histogram in the transverse momentum distribution 
represents the NLO pQCD predictions [S] listed in [15]. The experimental data are from 
CMS [T], ATLAS [2], DO [T5l[l6l|26] and CDF [I71|27]. 
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Figure 2: Different contributions to the total cross sections of top quark pair production in 
pp and pp collisions at the Tevatron and LHC energies. The dashed and dash-dotted curves 
correspond to the contributions from q*q* — )■ ti and g*g* — )■ ti subprocesses, respectively. 
The solid curve represents the sum of these components. The dotted curves correspond to 
the contribution of the sea quarks coming from the earlier gluon splittings. We use the 
KMR parton densities for illustration. The experimental data are from CMS [Ij, ATLAS [2], 
DO [T5l[l6l|26] and CDF [ni|27]. 



